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Decomposing Tensors into Frames 

Luke Oeding, Elina Robeva and Bernd Sturmfels 


Abstract 

A symmetric tensor of small rank decomposes into a configuration of only few vectors. 
We study the variety of tensors for which this configuration is a unit norm tight frame. 


1 Introduction 


A fundamental problem in computational algebraic geometry, with a wide range of appli¬ 
cations, is the low rank decomposition of symmetric tensors; see e.g. [D S Uni 1201 [2T]. If 
T = {tni 2 -id) Is a symmetric tensor in Sym^^lC”'), then such a decomposition takes the form 

r 

i=i 

Here Xj G C and = {vij,V 2 j, ■ ■ ■, Vnj) G C" for j = 1, 2,..., r. The smallest r for which a 
representation ([T]) exists is the rank of T. In particular, each is a tensor of rank 1. 

An equivalent way to represent a symmetric tensor T is as the homogeneous polynomial 


T 


E 


t 




Xi-, Xi. 


‘X. 
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If d = 2, then ([2]) is the identihcation of symmetric matrices with quadratic forms. Written 
as a polynomial, the right hand side of ([T]) is a linear combination of powers of linear forms: 

r 

T = '^Xj{vijXi+V2jX2-^ - \-VnjXnY. (3) 

i=i 

The decomposition in ([1]) and ([3]) is called Waring decomposition. When d = 2, it corresponds 
to orthogonal diagonalization of symmetric matrices. We could subsume the constants A, 
into the vectors Vj but we prefer to leave ([I]) and ([3]) as is, for reasons to be seen shortly. The 
(projective) variety of all such symmetric tensors is the r-th secant variety of the Veronese 
variety. The vast literature on the geometry and equations of this variety (cf. [18]) forms 
the mathematical foundation for low rank decomposition algorithms for symmetric tensors. 

In many situations one places further restrictions on the summands in ([T]) and (l3|), such 
as being real and nonnegative. Applications to machine learning in [1] concern the case 
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when r = n and the vectors vi,..., v„ form an ortho normal basis of M”. The article [21] 
characterizes the odeco variety of all tensors that admit snch an orthogonal decomposition. 
The present paper takes this one step fnrther by connecting tensors to frame theory 
We examine the scenario when the vj form a finite unit norm tight frame 
(or fnntf) of M"", an object of recent interest at the interface of applied fnnctional analysis 
and algebraic geometry. Consider a conhgnration V = (vi,...,Vr) G (M”)'' of r labeled 
vectors in R"". We also regard this as an n x r-matrix V = {vij). We call V a funtf if 

n 

V ■ = — ■ Id„ and vf, = 1 for i = 1, 2,..., r. (4) 

n 

i=i 

This is an inhomogeneons system of + r qnadratic eqnations in nr nnknowns. The funtf 
variety, denoted as in [B], is the snbvariety of complex affine space dehned by (| 1 |). 
For the state of the art we refer to the article [ 6 ] by Cahill, Mixon and Strawn, and the 
references therein. A detailed review, with some new perspectives, will be given in Section [ 2 l 
We homogenize the fnntf variety by attaching a scalar Aj to each vector Vj. The resnlt 
maps into the projective space P(Sym^(C"')) = of symmetric tensors, via the 

formnlas ([ID and ([3D- Onr aim is to stndy the closnre of the image of that map. This is 
denoted Tr,n,d- We call it the variety of frame decomposable tensors, or the fradeco variety. 
Here r, n, d are positive integers with r > n. For r = n, Tn,n,d is the odeco variety of |21j . 

Example 1.1. Let n = 3,(i = 4, and consider the symmetric 3x3x3x3-tensor 

T = b9{x\ + X2+x\) — 1 Q[x\x2 + X 1 X 2 + x\x^ + x\x2, + Xixl + X2x\) , . 

+ Q6{xlxl + xfxl + X^xl) + 96 (xIx2X3 -|- XiX^X^ -|- XiX2xl). ^ 


This ternary qnartic lies in 71 , 3 , 4 , i.e. this tensor has fradeco rank r = 4. To see this, note that 

T = ^{-5XI+X2+X3Y + ^{XI-5X2+X3)^+^{XI+X2-5X3)'^+^{3XI + 3X2 + 3X3)'^. ( 6 ) 
The corresponding fonr vectors, appropriately scaled, form a hnite nnit norm tight frame: 


V 




^ 7^4,3- 


(7) 


The fradeco variety 71 , 3,4 is a projective variety of dimension 6 and degree 74 in It is 
parametrized by applying rotation matrices p G SO 3 to all ternary qnartics of the form 

T = \i{-5xi + X2 + X3y + X2{xi-5x2 + X3y + \3{xi+X2-5x3)'^ + X4{3xi + 3x2 + 3x3Y. ( 8 ) 

The title of onr paper refers to the task of hnding the ontpnt (ED from the inpnt (ED. In this 
particular case, the decomposition can be found easily using Sylvester’s classical Catalecticant 
Algorithm, as explained in [20] Section 2.2]. In general, this will be more difficult to do. <) 
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The fradeco rank of a symmetric tensor T G Sym^(M"') is defined as the smallest r snch 
that T G Tr^n,d- This property does not imply that T also has a frame decomposition ([T]) 
of length r + 1. Indeed, we often have Tr,n,d 't- 'Tr+i,n,d- For instance, the odeco qnartic 
xf + X 2 + x\ lies in 73,3,4\74,3,4, by the constraint in Example 14.51 See also Example 13.51 
This paper is organized as follows. In Section [2] we give an introdnction to the algebraic 
geometry of the fnntf variety J>,n- This lays the fonndation for the snbseqnent stndy of 
fradeco tensors. Section [3] is concerned with the case of symmetric 2x2x • • • x2-tensors T. 
These correspond to binary forms (n = 2). We characterize frame decomposable tensors in 
terms of rank conditions on matrices. In Section 0] we investigate the general case n > 3, 
and we present what we know abont the fradeco varieties %,n,d- Section [5] is devoted to 
nnmerical algorithms for stndying Tr,n,d and for decomposing its elements into frames. 

2 Finite unit norm tight frames 

In this section we discnss varions representations of the fnntf variety J>,n- This may serve 
as an invitation to the emerging interaction between algebraic geometry and frame theory. 

Each variety stndied in this paper is dehned over the real held M and is the Zariski closnre 
of its set of real points. This Zariski closnre lives in affine or projective space over C. For 
instance, SO„ is the gronp of nxn rotation matrices p, and snch matrices have entries in R. 
However, when referring to SO„ as an algebraic variety we mean the irredncible snbvariety 
of dehned by the polynomial eqnations p ■ = Id„ and det(p) = 1. Likewise, a fnntf 

E is a real n x r matrix, bnt the fnntf variety lives in It consists of all complex 

solntions to the qnadratic eqnations (jl]). In the frame theory literatnre [5l [6], [121126] there is 
also a complex Hermitian version of J>,n, bnt it will not be considered in this paper. 

It is important to distingnish J>,n from the variety of Parseval frames, here denoted Vr,n- 
The latter is mnch easier than the former. The variety Vr,n is dehned by the matrix eqnation 

V-V^= Idn. 

The real points on Vr,n are smooth and Zariski dense, and they form the Stiefel manifold of 
all orthogonal projections R'’ —)■ R”. Hence Vr,n is irredncible of dimension nr — 

One featnre that distingnishes Vr,n from is the existence of a canonical map Vr,n+i 
Vr,n- Indeed, by Naimark’s Theorem [3], every Parseval frame is the orthogonal projection 
of an orthonormal basis of R'’, so we can add a row to H G Vr,n and get a matrix in Vr,n+i- 
There is no analogons statement for the variety We begin with the following resnlt. 

Theorem 2.1. The dimension of the funtf variety J^r,n is 

Ti 

dim(J>,n) = ~ 1) • ~ 1) provided r > n > 2. (9) 

It is irreducible when r > n + 2 > A. 

Proof. Cahill, Mixon and Strawn [6l Theorem 1.4] proved that is irredncible when 
r > n -|- 2 > 4. The dimension formnla comes from two articles: one by Dykema and Strawn 
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r 

n 

dim J>,n 

deg Rr,n 

# components & degrees 

3 

2 

1 

8-2 

8 components, each degree 2 

4 

2 

2 

12-4 

12 components, each degree 4 

5 

2 

3 

112 

irreducible 

6 

2 

4 

240 

irreducible 

7 

2 

5 

496 

irreducible 

4 

3 

3 

16-8 

16 components, each degree 8 

5 

3 

5 

1024 

irreducible 

6 

3 

7 

2048 

irreducible 

7 

3 

9 

4096 

irreducible 

5 

4 

6 

32-40 

32 components, each degree 40 

6 

4 

9 

20800 

irreducible 

7 

4 

12 

65536 

irreducible 


Table 1; Dimension and degree of the funtf variety in some small cases 


[121 Theorem 4.3(ii)] regarding the case when r and n are relatively prime, and one by Strawn 
[26l Corollary 3.5] which studies the local geometry for all r, n. In these articles it is shown 
that the real points in have a dense open subset that forms a manifold of dimension 
(n — 1) • (r — I — 1). The arguments in [6] show that the real points are Zariski dense in the 
complex variety Hence ([9]) is the correct formula for the dimension of J^r,n- HI 

Next to the dimension, the most important invariant of an algebraic variety is its degree. 
By this we mean the degree of its projective closure [HI §8.4], This can be computed using 
symbolic software for Grobner bases, or using numerical algebraic geometry software. The 
dimension and degree of for small r,n in Tabled] were computed using Bertini [2]. 

The case r = n + 1 is special. Here, the funtf variety decomposes into 2”+^ irreducible 
components, each of which is affinely isomorphic to the (”)"dimensional variety SO„. This 
will be explained in Corollary 12.101 The next example discusses one other exceptional case. 

Example 2.2 (r = 4, n = 2). Following ([4]), the dehning ideal of the funtf variety equals 
('^11 + '^12 + ~ 2, V 11 V 21 + V 12 V 22 + Vi‘j,V2i, + U 14 U 24 ) + (Di + ~ 1 f = 1, 2, 3,4). 

Note that this contains + V 22 + '*^23 + "*^24 “ 2. Using Grobner basis software, such as 
Macaulay 2 [T3], one checks that this ideal is the complete intersection of the six given 
quadrics, it is radical, and its degree is 48. Primary decomposition reveals that this ideal is 
the intersection of 12 prime ideals, each of degree 4. One of these associated primes is 

( — V22-, V12 + U21, U31 — U42, '^32 +'^41) + ( + G^2 ~ 1 • * = 1 ) 2 , 3 , 4 ). 

The irreducible variety of this particular prime ideal consists of the 2 x 4-matrices 

U = (i?i I R 2 ) , 
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where Ri and R 2 are rotation matrices of format 2x2. The other 11 components are obtained 
by replacing R^ by —Ri and permuting columns. The image of V under the map to binary 
forms is a linear combination of two odeco forms, one given by Ri and the other by i? 2 - ^ 

The real points of live in where = {m e M” : = 1} denotes the 

unit sphere. However, the vectors on these spheres will get scaled by the multipliers in 
([3]) when we pass to the fradeco variety %^n,d- To achieve better geometric properties and 
computational speed, we map each real sphere to complex projective (n—l)-space 

The projective funtf variety Qr,n is the image of Rr,n in To describe its equations, 

we use an n X r-matrix V = (vij) of unknowns as before, but now the i-th column of V 
represents coordinates on the i-th factor of We introduce the r x r diagonal matrix 

( n n n \ 

5^4,5^4, •••, ( 10 ) 

i=l i=l i=l ' 


The variety Qr^n is dehned by the following matrix equation: 

V ■ D-^ ( 11 ) 

n 

Each entry on the left hand side is a homogeneous rational function of degree 0. In fact, 
these functions are multihomogeneous: they dehne rational functions on (P"’“^)^. 

The challenge is to clear denominators in flTTl) . so as to obtain a system of polynomial 
equations that dehnes Qr,n as a subvariety of (P"'”^)'". Next we solve this problem for n = 2. 
For planar frames, equation flTT]) translates into the vanishing of the two rational functions 








— r 


and 






( 12 ) 


Consider the numerator of the rational function 


P-iQ = 


E 




2ivijV2j 


" 2 j 




Vij - V2jt 




E 


where i = \/ —1. 


Let P and Q denote the real part and the imaginary part of that numerator. These are two 
multilinear polynomials of degree r with integer coefficients in nii,ni 2 , • • • ,V 2 r- They dehne 
a complete intersection, and, by construction, this is precisely our funtf variety in (P^)^: 

Lemma 2.3. The projective funtf variety Qr ,2 is a complete intersection of codimension 2 in 
(P^)^, namely, it is the zero set of the two multilinear forms P and Q. 

Here are explicit formulas for the multilinear forms that dehne Qr ,2 when r < 5: 

Example 2.4. If r = 3, then P = 3 vuVi 2 Vi 3 + VUV 22 V 23 + V 21 V 12 V 23 + V 21 V 22 V 13 and Q = 
vnVi2V23+viiV22Vi3+V2iVi2Vi3+3v2iV22V23- If r=4, then P = 4:{vuVi2Vi3Vu-V2iV22V23V24:) and 

Q = 2viiVi2Vi3V2i + 2niini2n23t'i4 + 2 ^ 11 ^ 22 ^ 13^14 + 2t;iin22n23^’24+ 

2n2l'*^12'i^l3'yi4 + 2t;2l'?^12'i^23'*^24 + 2n2l'?^22'i^l3'i^24 + 2n2l'i^22'*^23'i^l4- 
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If r = 5, then 


P = ^ViiVi2VlzVlAVi5 - ViiVi2Vi‘iV24V2b - ViiVi2V2^Vi4V25 - ViiVi2V21,V24Vl5 

-ViiV22Vl‘iViiV2^ - ViiV22Vi:iV2AVl5 - ViiV22V2‘iVlAVi^ - ?>ViiV22V2‘iV2iV2h 
-V2 iVi2Vi^ViaV2^ - ^21 Wl2t'l3^^24Wl5 ” W2lt'l2W23^'l4Wl5 “ 3n2lt'l2W23^'24^^25 
—V2lV22Vl‘iViAVi^ — '3tV2lV22Vl'iV2AV2b ~ 3n2in22'y23'i'l4'y25 ~ 3n2in22'y23'i'24'yi5, 


and Q is obtained from P by switching the two rows oiV. <) 

Such formulas are useful for parametrizing frames. We write the equations for Qr ,2 as 


P\ ^ fmii mi 2 \ /uiA ^ 

Q 1 \^m 2 i m22 ) V^ 2 r/ W / 


The matrix entries mij are multilinear forms in (un : n 2 i), (^^12 : '^ 22 ), • • •, (wi,r-i : V 2 ,r-i)- 
Using the quadratic formula, we solve the following equation for one of its unknowns: 


"iiim 22 = mi2m2i. 


(13) 


This dehnes a hypersurface in from which we can now easily sample points. The 

point in the remaining rth factor P^ is then recovered by setting vir = mi 2 , V 2 r = —''Tin. 

For n > 3, we do not know the generators of the multihomogeneous prime ideal of Qr,n- 
Here are two instances where Macaulay2 [13] succeeded in computing these ideals: 

Example 2.5. The variety Qa,z is a threefold in (P^)^. Its ideal is generated by 34 quartics. 
Among them are the equations that dehne the six coordinate projections into (P^)^, like 


8 (^iWi 2 + ^ 2^22 + ^ 3 ^ 32 ) + 18(nlln21^^12^^22 + ViiV‘iiVi2V‘i2 + V2lV‘iiV22V‘i2) 




- UTi V 


11 *^32 


- 


- V2lVh - V 


31 


^^ 1^2 


— V 


31 


^' 22 - 


Example 2.6. Let r = 5 and n = 3. By saturating the denominators in (fTT]) . we found 
that the ideal of ^ 5^3 is generated by a 120-dimensional SOa-invariant space of sextics. The 
following polynomial (with 60 terms of Z^-degree ( 2 , 2 , 2 , 0 , 0 )) is a highest weight vector: 




4 - 




+iStijj^ 1122'y3211231133 


23 + 5«ii« 1243 + 45i;hlll2l’22l'13l^23 

5’^1i42'*'13 ~ ^’^1i4243 + ^^ll'^32'^33 
+45))llD21'yi211221l23 + 18i;ilD21'W12'«22*l33 + 27 dii1>21' 1 ' 12 1132'y23'1133 
+45i;iltl3ltlf2'!113'!l33 + 27uiltl3lUl2ll22l’23l'33 + 45lllli;3li;i2D32'l'f3 
+ 18llllD3lll|2l’13l'33 + 27i;ili;3li;22'l’32l’13ll23 + 45'illli;3ltl32l>13l’33 


45'ilfji;i2t>32l’13'!l33 
45DllD211'f2'''131123 


■ 5l'lll'22l'l3 ■ 


- 5lljjll22l’23 4:llhll22l’33 


451111112111121’22 Ilf 3 + 18llll1i211l32l’131'23 
'22III3II23 + 27llll1l211’221'321’131’33 
18llllll3llll2ll32l’|3 + 451111 II 31 II 12 II 32 II 33 - 4ll2il’f2''33 
^3 + 5ll|illf2«|3 + 45ll|j_1Il2ll22lll3_l’23_ 


SllilllL" 


+ 18ll^llll2ll32lll3l'33 + SlliillljllL + 50414243 + ^'>’21^h^33 + 45ll|ill22ll32ll23l'33 “ 4llii42^13 + ^'>’21^32^M + 

+ 18iI21113i 42''231’33 + 27lI211'311’12112211131’33 + 27ll211’311’121'321'131’23 + 45ll211l311'22''231’33 + 18ll211’311’221’321’i3 
+45ll211’311i221'321l23 + 4:5V21V33V22V32V^^ + 45li211’3l42''231133 + ^'>^ 31 ^ 12^13 “ ‘i^31^12^^3 + 54l''l243 + 184l'''121i221’131123 

5ll?, 1l2o4i + 5O1I0, 11?o1l?,. 


■ 54i''i2''i3 ^4i''i2''J3 

+45ll§j^ 1112113211131133 - 44 i4243 + ^^31^22^23 + 54i 4ollh + 4511?, 1122113211231133 + 5ll. 


■^31 “^22 *^33 


“31 “32 “13 


The ideal of ^ 5^3 C (P^)^ has 10 generators like this, each spanning a one-dimensional graded 
component. It has 30 components of degrees like (2, 2,1,1, 0), each generated by a polynomial 
with 78 terms. Finally, it has hve 16-dimensional components of degrees like (2,1,1,1,1). <0 
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In order to sample points from the funtf variety we can also use the following 

parametrization found in [3122]. We write V = {U', W), where U' is an n x n-matrix and W 
is an {r — n) X n-matrix. For the columns of W we take arbitrary points on the unit sphere 
In practice, it is convenient to hx a rational parametrization of so as to ensure 
that W has rational entries Wij. For instance, for n = 3 we might use the following formulas: 


Wij 


2Aj/ij 


W2j 


2XjUj 




W3j 






2 ’ 
3 


where G 1>. (14) 


After these choices have been made, we £x the following n x n-matrix with entries in 

r 


S = 


■Id. 


W-W^ 


n 


(15) 


It now remains to study all n x n-matrices 1/ = (uij) that satisfy 


[/■D ^ -[/^ = S, where D = diag( , ..., 

2=1 2 = 1 


For any such [/ we get a funtf V = {U\ W) G by setting U' = U ■ D For random 
choices in (iMj), the matrix S is invertible, and the previous equation is equivalent to 

D = U^-S-^-U. (16) 

This identity of symmetric matrices dehnes ("' 2 ^) equations in the entries Uij of U. The 
equation in position (i, j) is bilinear in {uu, U 2 i, ■ ■ ■, u^i) and (uy, U 2 j -,..., Unj)- We solve the 
system ffT6l) iteratively for the columns of U. We begin with the (1,1) entry of ffTHll . There 
are n — 1 degrees of freedom to £11 in the £rst column of U, then n — 2 degrees of freedom to 
£11 in the second column, etc. This involves repeatedly solving quadratic equations in one 
variable, so each solution lives in a tower of quadratic extensions over Q. In summary: 


Proposition 2.7. 


The equations jlly , 11.5]/ . ( lih)) represent a parametrization 


'-'j ^ r.Ti • 


The rotation group SO„ acts by left multiplication on the funtf variety J^r,n- There is a 
natural way to construct the quotient as an algebraic variety, namely by mapping 

it into the Grassmannian Gr(? 7 ,, r) of n-dimensional subspaces of C^. This is described by 
Gahill and Strawn in O Section 3.1], and we briehy develop some basic algebraic properties. 

We here dehne Gr(n, r) to be the image of the Pliicker map —?■ that takes 

an n X r-matrix V to its vector p = p{V) of n x n-minors. The coordinates pi of p are 
indexed by the set of n-element subsets of [r] = {l,2,...,r}. With this dehnition, 

Gr(n, r) is the affine subvariety of dehned by the quadratic Pliicker relations, such as 
P12P34 — P13P24 + P14P23 = 0forn = 2,r = 4 . The dimension of Gr(n, r) is (r — n)n + 1. Note 
that if VV'^ = (r/n) ■ Id„, then the Cauchy-Binet formula (cf. 0 Prop. 6 ]) implies 


= 



( 17 ) 


The real points in Gr(? 7 ,, r), up to scaling, correspond to n-dimensional subspaces of 
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Proposition 2.8. The image of J^r,n under the Pliicker map is an affine variety of dimension 
(r—n)n — r + 2 in the Grassmannian Gr(?7,, r) C . It is defined by the equations 

= (-) for i = 1 , 2 ,..., r. (18) 

I:iel 

The real points in this image correspond to SOn-orbits of n-dimensional frames in J>,n- 

Note that adding up the r relations in ffTSj) and dividing by n gives precisely flT7|) . 

Proof. Both lFr,n and the constraints flTSl) are invariant under SO„. Suppose that V G 
satishes VV'^ = {r/n) ■ Id„. We may assume (modulo SO„) that the i-th column of V is 
(a, 0,, 0)^ for some a G C. Let V be the matrix obtained from V by deleting the hrst 
row and i-th column. Then V ■ = {r/n) ■ Idn_i. Any pi with i G / equals a times the 

maximal minor of V indexed by Applying (lT7j) to V, this gives 



I:iel 


Hence (1T8|) holds if and only if a = ±1, and this holds for all i if and only if V lies in J>,n- 
The dimension formula follows from Theorem 12. II because SO^ acts faithfully on J>,n- CH 

Example 2.9. Let n = 2. If r = 5, then our construction realizes ^^ 5 , 2/802 as an irreducible 
surface of degree 80 in Its prime ideal is generated by the ten quadratic polynomials 

PIAP23 - P13P2A + P12P34, P15P23 - P13P25 + P12P35, P15P2A - PIAP25 + P12P45, P 15 P 3 A - P1AP35 
+P13P45, P25P34 - P24P35 + P 23 PA 5 , Pu + PI3 + Pl 4 + Pl 5 - 5 / 2 , Pu + pIs + pIa + Plb “ 5 / 2 , 
P ?3 + PI3 + pIa + P35 - 5 / 2 , pIa + pIa + pIa + P% - 5 / 2 , pl^ + pl^ + pl^ + p% - 5 / 2 . 

If r = 4, then ^^ 4 , 2/802 is a reducible curve of degree 24 in C®. Its dehning equations are 

P 14 P 23 -P 13 P 24 +P 12 P 34 = 0, Pi2+P?3+Pi4 = PI2+PI3+P2A = P13+P23+P3A = pIa+pIa+pIa = 2- 

As in Example 12.21 this curve breaks into 12 components. One of these 12 irreducible curves 
is { p G C® : P 12 = P 34 = 1, Pl 3 = P 2 A, PlA = -P23, pIs + ^24 = 1} • 0 

The analogous decomposition is found easily for the case r = n + 1. Here, there are no 
Pliicker relations, so Gr(n, n + 1) ~ For convenience of notation, we set g* = p[n+i]\{j} in 
(IT^ . The quotient space J^n+i,n/SOn is the subvariety of dehned by the equations 

?? + ^2 + ■ • • + + qI+i = {n + + <li for i = 1 , 2 ,..., n + 1 . 

These are equivalent to the following equations, which imply Gorollary 12.101 

q! = ql = ql = --- = qUi = + 

Corollary 2.10. The quotient space J>i+i,n/80„ is a variety consisting of 2”"''^ isolated 
points in = Gr(n, n + 1), namely those points with coordinates ±(n+ 






Any of the components of Tn+\,n can be used to parametrize our variety Tn+\,n,d- 

Example 2.11. Let n = 3. The point p = \/3(|, |, |, |) in Gr(3,4) corresponds to the 
SOa-orbit of the frame V in Example II. 1 [ The variety ^ 4^3 can be parametrized as follows: 




1 — 2 y‘^ — 2z‘^ 2xy — 2zw 2xz + 2yw 


■3 1 1 -5' 


2 xy + 2zw 1 — 2 x‘^ — 2z^ 2yz — 2xw 


3 1-51 


2 xz — 2yw 2yz + 2xw 1 — 2 x‘^ — 2y^ 


3-511 



0 0 0 

z/2 0 0 

0 1/3 0 

0 0 U4 


The 3 X 3-matrix on the left is the familiar parametrization of SO 3 via unit quaternions. 
This gives the parametrization of the fradeco variety 71, 3 , d seen in ([8]). <) 

The embedding of 7>,n/SOn into Gr(r, n) via (flSD connects frame theory with matroid 
theory. The matroid of V is given by the set of Pliicker coordinates pi that are zero. If all 
Plucker coordinates are nonzero, then the matroid is uniform. It is a natural to ask which 
matroids are realizable over M when the additional constraints flTSll are imposed. 

The discussion in [5], Section 3.2] relates frame theory to the study of orhitopes [22]. Gahill 
and Strawn set up an optimization problem for computing Parseval frames that are most 
uniform. Their formulation in [5l p. 24] is a linear program over the Grassmann orbitope, 
which is the convex hull of Gr(n,r) intersected with flTTl) . The same optimization problem 
makes sense with Gr(n,r) replaced by 7>,„/SO„, or, algebraically, with f[T7|) replaced by 
ffTS]) . If n = 2 , then the former problem is a semidefinite program. This is the content of \T2l 
Theorem 7.3]. For n > 3, the situation is more complicated, but the considerable body of 
results coming from calibrated manifolds, such as [22l Theorem 7.5], should still be helpful. 


3 Binary forms 

We now commence our study of the fradeco variety %^n,d- In this section we focus on the 
case n = 2 of binary forms that are decomposable into small frames. The case r = 2 is the 
odeco surface known from [211 §3]- Proposition 3.6 in [21] gives an explicit list of quadrics 
that forms a Grobner basis for the prime ideal of 12 , 24 -, these are here expressed as the 
2 X 2-minors of a certain 3 x (d —3)-matrix M. 4 . What follows is our main result in Section [31 
We are using coordinates (to : ■ ■ ■ : td) for the space = P(Sym^(C^)) of binary forms. In 
the notation of ([2]), the coordinate ti would be tiii...i 222 - -2 with i indices 1 and d — i indices 2. 

Theorem 3.1. Fix r G {3,4,..., 9} and d > 2r — 2. There exists a matrix Mr such that: 

(a) Its maximal minors form a Grobner basis for the prime ideal of %, 24 - 

(b) It has r — 1 rows and d — r + 1 columns, and the entries are linear forms in to,... ,td. 

(c) Each column involves r of the unknowns U, and they are identical up to index shifts. 

These matrices can be chosen as follows: 

_ fto — 3^2 ti — 3to t2 — 3^4 G — 3^5 • ■ ■ td-3 — 3td-l 

^ \3ti — G 3^2 — G SG — ts 3G — to ■ ■ ■ 3G-2 ~ G 
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/ to + ^4 

tl P tg 

t2 + tg 

tgptj ■ ■ ■ 

td—i T tj;; \ 


AiA = 

ti-tg 

t2 — ti 

ts — tg 

ti — tg ■ ■ ■ 

td-3 + td-1 1 

(20) 


\ t2 

ts 

ti 

tg 

td-2 ) 



^tg p 5^2 

tl P 5tg 

^2 + 5^4 

ts P 7)tg ■ ■ ■ 

td-g + 5t(i_3^ 


Ai5 = 

ti - Stg 
St2 - t4 

t2 — Sti 
3^3 ~ tg 

ts ~ 3^5 

3^4 — tg 

ti — stg ■ ■ ■ 
Stg — tj ■ ■ ■ 

td—A 3t^_2 

3td-3 — td_l 

(21) 


^5^3 -|- tg 

5ti p tg 

5tg p tj 

Stg -|- tg • • • 

5td_2 P td / 



^tg p 3^2 

tl P Stg 

^2 + 3^4 

tg P Stg ■ ■ ■ 

td-g + iitd-i\ 



tl + tg 

t2 + tg 

tg -|- tj 

tiptg 

td-g + td_i 


Aie = 

t2 — ti 

ts — tg 

ti — tg 

tg — tj 

td-A — td -2 

(22) 


ts 

U 

tg 

tg 

td -3 



^3^4 -|- tg 

Stg -|- ^7 

3^6 + tg 

Stj -|- tg • • ■ 

3trf_2 PtdJ 


The first column of Ai 7 is {Stg -|- 7 ^2, ti + 5tg, t 2 — Si 

4 , 3t3 — tg, 5ti 

+ tg, 7tg P StjY', the first 


column of Aig is (to + 2 ^ 2 , ti + Sta, tg — tg, t 2 + te, Sts + tr, ‘^te + ts)^, the first column 
of Aig is (5to + 9t2, 3ti + 7tg, t 2 + 5t4, tg — Sts, 3t4 — tg, 5ts + tj, 7tg + Stg, Qty + 5tg)'^. 

We conjecture that the same result holds for all r, and we explain what we currently 
know after the proof. Let us begin with a lemma concerning the dimension of our variety. 


Lemma 3.2. The fradeco variety Tr^ 2 ,d is irreducible and has dimension min(2r — 3,d). 

Proof. For d> 5, the funtf variety J >,2 C (S^)^ is irreducible, by Theorem 12.11 and hence so 
is its closure Qr ,2 in (P^)^. While the two special varieties ^^ 3,2 and Jvi ,2 are reducible, the 
analyses in Example 12.21 and Corollary 12.101 show that ^ 3^2 and ^ 4,2 are irreducible. 

Regarding Qr ,2 as an affine variety in we obtain 7^,2,d as its image under the map 


ti = + VI 2 V 22 ' 


+ '^13^23 




for i = 0,1, ■ ■ ■ ,d. 


(23) 


This proves that Tr^ 2 ,d is irreducible. To see that it has the expected dimension, consider 
the r-th secant variety of the rational normal curve in P'^, which is the image of the map 
([;; 2 xr pd gjyg]^ py ([23]). It is known that this secant variety has the expected dimension, 
namely min(2r — l,d), and the hber dimension of the map (|2^ does not jump unless some 
2 X 2-minor of R = (vij) is zero. Since codim(^j,, 2 ) = 2, by Lemma iTSl the claim follows. □ 


Proof of Theorem \3.1[ We hrst show that the maximal minors of our matrices Air vanish 
on the fradeco variety Tr^ 2 ,d for r = 3,4,..., 9. After substituting the parametrization (l23l) 
for tg, ^ 1 , • • •, td, we can decompose these matrices as follows: 


where 


Air = Mr 


hP 

n.d-r-l 

Al ^21 

^.d—r 

. d.—r—l . 

.d-r-2 .2 

V 12 

^12 

V 22 

A 2 

^22 

^.d—r 

d—r — \ 

d—r—2 2 

\filr 

Vir 

V2r 

Vir 

V2r 


vir\ 

vir 


d-r 

hr 




({VI 2 - 3vlfi)V2l {VI 2 - 3vl2)v22 - 3uf3)u23\ 

V(3u|i - vlfijVii {3vl2 - vl 2 )vi 2 (3u|3 - vlg)vig) 
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Ms 


/ V 21 + v\i V 22 + '^12 '^23 + "^24 + \ 

M 4 = - WnV 21 t'l 2 wf 2 - ^ 12^22 ^13^23 “ ^13^23 t'l4t’i4 “ ^14^24 , 

V VI 2 VI 2 vl,vl, vl,vl, ) 


( 


^21 


+ 5u? 


^22 


+ bvl 


^23 


+ 54 


13 

,3 


^^24 

,,4 


5^4 


■^25 


+ 545 \ 

.,3 " 


Ull^l - V 12 VI 2 - 31112^22 ^13'(^23 “ 3lll3^^23 ''^1444 “ 3^1444 ^^15^^25 “ 


3^2 „3 


U 91 - 


■^ 11^^21 


V 54l4l + V 


'^111121 


3^2 „3 


VrfO — 


^ 12^^22 


fl24'22 


3^^13^^23 - 


^11 


54242 + ^ 


12 


5^^13^^23 + '^' 


V 13 V 23 
5 


3,,2 „.3 
0(114 (124 


13 


511441124 + V 


VU'<^2A 
5 


14 


3^^1545 - ^^15^25 

54545 + 45 / 


and similarly for Mg, My, Mg and Mg. We claim that the matrices M^ have rank < r — 1 
whenever V G J>, 2 - Equivalently, the (r — 1) x (r — 1) minors of M^ lie in the ideal of Qr, 2 - It 
suffices to consider the leftmost such minor since all minors are equivalent under permuting 
the columns of V. For each r < 9, we check that the determinant of that minor factors as 


(miim 22 - mi 2 m 2 i) • JJ 4ii'y2j -'^ 2 *'^^), (24) 

l< 2 <j<r —1 

where the left factor is the polynomial of degree 2r — 2 given in flTdl) . That polynomial 
vanishes on Qr, 2 - This implies rank(Mj,) < r—2 on Qr, 2 , and hence rank(A4r) < r—2 on Tr^ 2 ,d- 
Fix the lexicographic term order on C[fo, 4,..., fj. We can check that, for each r G 
{3,4,..., 9}, the leading monomial of the leftmost maximal minor of Air equals fo^ 2^4 • • • tr- 2 - 
Hence all Kiaximal minors of Air are squarefree, and they generate the ideal 

Ir^d ^ ^21 ^( 2^(3 * * * i^r—l * 2 + ^ 1+2 + ^ 2 : ^ 2+2 + ^3; ^3 +2 + ^ 4 , • • • , tr—2A 2 + ^r —1 + d 2 

This squarefree monomial ideal is pure of codimension d — 2r + 3 and it has degree (^^4^• 
This follows from [191 Theorem 1.6]. Indeed, in Murai’s theory, our ideal A^d is obtained from 
the power of the maximal ideal by applying the stable operator given by a = ( 2 ,4, 6 ,...). 
Combinatorial analysis reveals that the ideal Ir,d is the intersection of the prime ideals 

{^jo y ^ji y ^32 y^jiy ■ ■ ■ y ^jd- 2 r +2 ) > 

where Jo,4,24,... are even, ji, js, J 5 ,... are odd, and 0 < jo<ji<j 2 < • • • <jd- 2 r +2 < d. Note 
that number of such sequences is (^J^r+s) ~ • H^^^e the codimension and degree of 

Ir,d are as expected for the ideal of maximal minors of an (r —1) x (d—r+l)-matrix with linear 
entries m Ex. 19.10]. The monomial ideal Ir,d is Cohen-Macaulay because its corresponding 
simplicial complex is shellable (cf. [2ll §111.2]). Indeed, if we list the associated primes in a 
dictionary order for all sequences joji 22 • • • jd- 2 r +2 as above, then this gives a shelling order. 

Using Buchberger’s S-pair criterion, we check that the maximal minors of Air form a 
Grobner basis. We only need to consider pairs of minors whose leading terms share variables. 
Up to symmetry, there are only few such pairs, so this is an easy check for each fixed r < 9. 

Since A^d is radical of codimension d — 2r -|- 3, we conclude that the ideal of maximal 
minors of Air is radical and has the same codimension. However, that ideal of minors is 
contained in the prime ideal of A, 2 ,d, which has codimension d — 2r -|- 3 by Lemma [3.21 
Therefore, we now know that %, 2 ,d is one of the irreducible components of the variety of 
maximal minors of Air- To conclude the proof we need to show that the latter variety is 
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irreducible, so they are equal. To see this, we £x r and we proceed by induction on d. For 
d = 2r — 2, when Mj. is a square matrix, this can be checked directly. To pass from d to 
d + 1, we factor the matrix as Mj. times the rank r Hankel matrix associated with a funtf V. 
Increasing the value of d to d + 1 multiplies the i-th row of the Hankel matrix by vu and it 
adds one more column. This gives us the value for the new variable t^+i. Now, since that 
variable occurs linearly in the maximal minors, its value is unique. This implies that the 
unique rank r — 2 extension from the old to the new Air must come from the funtf V. □ 

We established Theorem 13.11 assuming that r <9, but we believe that it holds for all r: 

Conjecture 3.3. For allr > 3 there exists a matrix Air which satisfies the properties (a),(h) 
and (c) in Theorem \3.1\ When r is odd, the matrix Air is given by 

/^2 ^3 • • • id-r+2\ 

*3 *4 • • • td_r+3 

\tr tr+1 • • • td / 

V r/ V r —4/ 

We do not know yet what the general formula for Air should be when r is even. The 
following systematic construction led to the matrices Mr and Air in all cases known to us. 
Let P and Q be the multilinear forms in Lemma 12.31 that dehne Qr, 2 - Let Fj denote the 
polynomial of degree 2r — 2 obtained by eliminating vij and V 2 j from P and Q. Let Gj 
denote the product of all (^ 2 ^) minors VikV 2 i — ViiV 2 k of V where j ^ {k,l}. Each product 
FjGj is a polynomial of degree r(r — 1). Note that Fr is miim 22 — mi 2 m 2 i in flT^ . and 
FrGr is fl2T|) . Now, the ideal {FiGi, F 2 G 2 , ■ ■ ■ ,FrGr) is Cohen-Macaulay of codimension 2. 
By the Hilbert-Burch Theorem, the FjGj are the maximal minors of an (r — 1) x r-matrix 
Air, which can be extracted from the minimal free resolution of {FiGi ,..., FrGr)- This is 
precisely our matrix. In order to extend Theorem 13.11 and to End the desired Air for even r, 
we need that all entries of the Hilbert-Burch matrix Mr have the same degree r. 

Remark 3.4. The singular locus of 7^,2,d is dehned by the (r—2) x (r—2)-minors of Air- It 
would be interesting to study this subvariety of and how it relates to singularities of Pr, 2 - 
For instance, for r = 4, this singular locus is precisely the odeco variety 72, 2 ,d, and, using 
Proposition 3.6], we can see that its prime ideal is generated by the 2 x 2-minors of Ai^- 

In Section E] we shall see how the matrices Air can be used to End a frame decomposition 
of a given symmetric 2x2x • • • x2-tensor T. We close with an example that shows how this 
task difiers from the easier problem of constructing a rank r Waring decomposition of T. 

Example 3.5. Let r = 4, d = 8, and consider the sum of two odeco tensors 

T = + {x — y)^ + {x -\- y)^ = 3x^ + 56x^y^ + liQx'^y'^ + 56x^y^ + 3y^- 

The coordinates of this tensor are to = t^ = 3, t 2 = t^ = to = 2, and G = to = H = P = 0. 
Here, the 3 x 5-matrix Ai^, has rank 2. This veriEes that T lies in 74 , 2 , 8 ; in accordance with 
Example 12.21 However, the 4 x 4-matrix Ai^ is invertible. This means that T does not lie 
in 75,2,8- In other words, there is no funtf among the rank 5 Waring decompositions of T. <) 


/ ^0 


Mr-2 


*1 

t2 


\ 

td-r+l 


*d-2 / 
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4 Ternary Forms and Beyond 


We now move on to higher dimensions n > 3. Our object of study is the fradeco variety 


rr,n,d C P(Sym,(C")). 


A very hrst question is: What is the dimension of Tr^n,d ? In Lemma 13.21 we saw that 
dim(7j^ 2 ,d) = 2r — 3. The following proposition generalizes that formula to arbitrary m 

Proposition 4.1. For all r > n and d> 3, the dimension of% n,d is bounded above by 


min 


(n — l)(r — n) + 


[n 


l)(n 

2 


2 ) 


+ r - 1, 


hn + d — 1\ 

V d ) 



(26) 


Proof. The right number is the dimension of the ambient space, so this is an upper bound. 
The left number is the dimension of x by the formula in Theorem l2.ll The formula 
(Pj) expresses our variety as the (closure of the) image of a polynomial map 

J^r,n X P''"^ -)■ Tr,n,d- (26) 

The dimension of the image of this map is bounded above by the dimension of the domain. □ 


We conjecture that the true dimension always agrees with the expected dimension; 

Conjecture 4.2. The dimension of the variety %^n,d is equal to / f^) for allr > n and d >3. 

This conjecture is subtler than it may seem. Let arJydP"'~^ denote the Zariski closure 
of the set of tensors of rank < r in P(Sym^(C"')). Geometrically, this is the r-th secant 
variety of the d-th Veronese embedding of P"“^. It is known that o'rr'rfP””^ has the expected 
dimension in almost all cases. The Alexander-Hirschowitz Theorem (cf. [HIlH]) states that, 
assuming d > 3, the dimension of o'rr'dP””^ is lower than expected in precisely four cases: 

(r,n,d) e {(5,3,4), (7,5,3), (9,4,4), (14,5,4)}. (27) 

One might think that in these cases also the fradeco subvariety Tr,n,d has lower than expected 
dimension. However, the results summarized in Theorem |4]3] suggest that this is not the case. 

Theorem 4.3. Consider the fradeco varieties Fr,n,d in the cases when n > 3 and 1 < 
dim(Tr^n,d) ■ codim(7j^^„^rf) < 100. Table\B gives their degrees and some defining polynomials. 
The last column shows the minimal generators of lowest possible degrees in the ideal ofTr,n,d- 

Computational Proof. The dimensions are consistent with Conjecture 14.21 They were verihed 
by computing tangent spaces at a generic point using Bertini and Matlab. The degrees were 
computed with the monodromy loop method described in Subsection 15.2.11 The numerical 
Hilbert function method in Subsection 15. 2. 21 was used to determine how many polynomials of 
a given degree vanish on Tr^n,d- This was followed up with computations in exact arithmetic 
in Maple and Macaulay2. These conhrmed the earlier numerical results, and they enabled us 
to hnd the explicit polynomials in Q[T] that are listed in Examples 14.4114751 and WM In the 
cases where we report no equations occurring below a certain degree, this is a combination 
of Corollary 14.101 and the numerical Hilbert function computation. □ 
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variety 

dim 

CO dim 

degree 

known equations 

74,3,3 

6 

3 

17 

3 cubics, 6 quartics 

74,3,4 

6 

8 

74 

6 quadrics, 37 cubics 

74,3,5 

6 

14 

191 

27 quadrics, 104 cubics 

75,3,4 

9 

5 

210 

1 cubic, 6 quartics 

75,3,5 

9 

11 

1479 

20 cubics, 213 quartics 

76,3,4 

12 

2 

99 

none in degree < 5 

76,3,5 

12 

8 

4269 

one quartic 

77,3,5 

15 

5 

> 38541 

none in degree < 4 

78,3,5 

18 

2 

690 

none in degree < 5 

7io,3,6 

24 

3 

> 16252 

none in degree < 7 

75,4,3 

10 

9 

830 

none in degree < 4 

76,4,3 

14 

5 

1860 

none in degree < 3 

77,4,3 

18 

1 

194 

one in degree 194 


Table 2 : A census of small fradeco varieties 


We shall now discuss some of the cases appearing in Theorem 14.31 in more detail. 
Example 4.4. The 6-dimensional variety 74 , 3,3 C P® has the parametrization 

tsoo = + ^13 + ^14, 

^030 = '^21 + '^22 + '^2‘i + "^24) 

too3 = + VI 2 + nfs + n|4, 

^012 = V 2 W 3 I + ^22^32 + ^23^33 + ^24t'34) 

^021 = '*^21'*^31 + ^22^32 + + V 24 V 34 , 

tl02 = VlWil + Vi2vi2 + 

tl20 = V11V21 + V12V22 + V13V23 + ViiV24^, 

t201 = vllVsi + VI2V32 + vl^V33 + ^ 14 ^ 34 , 

^210 = vl^V21 + VI 2 V 22 + vl^V23 + Wm'^24, 

till = 'yil'y2l'y31 + V12V22V32 + 'yi3'?^23'i^33 + ni4'y24'i'34- 


Here the matrix V = {vij) is given by the parametrization of ^ 4,3 seen in ([ 8 ]) of Example ll.il 
Using exact linear algebra in Maple, we hnd that the ideal of 74 , 3,3 contains no quadrics, 
but it contains three linearly independent cubics and 36 quartics. One of the cubics is 

0*123 + 20*145 + 20*345 — 0*126 ~ 0*236 ~ 40*456, (29) 


where Cijk denotes the determinant of the 3x3 submatrix with columns i, j, k in 

( tsOO ^210 ti20 ^201 till ti02 

t210 tl20 to30 till to21 toi2 

t201 till to21 ti02 toi2 to03 

The other two cubics are obtained from this one by permuting the indices. The resulting 
three cubics dehne a complete intersection in P®. However, that complete intersection strictly 
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contains 74 , 3,3 because the three cubics have only 30 multiples in degree 4, whereas we know 
that 36 quartics vanish on 74 , 3 , 3 . Using Macaulay2, we identihed six minimal ideal generators 
in degree 4, and we found that the nine known generators generate a Cohen-Macaulay ideal 
of codimension 3 and degree 17. Using Bert ini, we independently verihed that fradeco 
variety 71 , 3,3 has degree 17. This implies that we have found the correct prime ideal. <C> 

Example 4.5. The variety 71 , 3,4 is also 6 -dimensional, and it lives in the of ternary 
quartics. The parametrization is as in (|28|) but with quartic monomials instead of cubic. 
Among the ideal generators for 71 , 3,4 are six quadrics and 37 cubics. One of the quadrics is 


8(^013 ^ 004 ^ 022 ) + 8(^031 ^ 022 ^ 040 ) + 8(^211 t 202 ^ 220 ) + I8{t 

+ (f004 ^O4O + 19fo22~20toi3^ 03l) + (^004f220+^022O02~2t013Oll 


112“Uo 3U21 ) -|- 18(t^21 ^ II 2 U 30 ) 
) -|- (to40^202+^022^220~2fo3lOll)- 


A Bert ini computation suggests that the known generators suffice to cut out 71 , 3 , 4 . We 
also note that the 27 quadrics for 71 , 3,5 come from the 6 quadrics for 71 , 3 , 4 . For instance, 
replacing each variable Ujk by yields the quadric 8^014 + 8^932 + ■ • ■ +19^023 for 71 , 3 , 5 . ^ 

Example 4.6. The fradeco variety 71 , 3,4 is especially interesting because (5,3,4) appears 
on the Alexander-Hirschowitz list (|27ll . The unique cubic that vanishes on 71 , 3,4 is 

46^022^202^220 + 73^112^121^211 — 4:^004^040^400 + 19 [^013^ 13of30l]2 — 50[fo04fll2 ]3 22 (^ 004 ^ 220 ] 3 

— 18[to22f 21 i] 3 + 50[fo04fo22f202]3 + 26[fo04fl3of31o]3 + 100[^013^ 103^ 112]3 ~ 53[^013^ 12lf31o]3 
+ 5 (^ 004 ^ 022 ^ 400 ]6 ~ 50(tQ;^3f202]6 ~ ^[toist22o]6 + 4:5(fo04fo3lf21l]6 ~ 4:0(fo22f202]6 + 5(to04fo22f22o]6 

+40(^022^ 112]6 ~ 5(fo04fl3o]6 ~ 45(to04fl2l]6 ~ 10(to04fll2fl3o]6~45(toi3fo22f21l]6+35(toi3fo3lf202]6 
+ 10(toi3fl03fl3o]6 + 10(toi3fll2fl2l]6 ~ 80(^013^ 112f30l]6 + 80(foi3f202f21l]6 + 8(foi3f21lf22o]6- 


This polynomial has 128 terms: each bracket denotes an orbit of monomials under the Ss- 
action, and the subscript is the orbit size. In addition, six fairly large quartics vanish on 71 , 3 , 4 . 
The seven known generators cut out a reducible variety of dimension 9 in P^^. The fradeco 
variety 71 , 3,4 is the unique top-dimensional component. But, using Bert ini, we found two 
extraneous components of dimension 7. Their degrees are 120 and 352 respectively. <0 

We close this section by examining the geometric interplay between fradeco varieties 
and secant varieties. We write ari2d^^~^ for the r-th secant variety of the d-th Veronese 
embedding of This lives in P(Sym^(C")) and comprises rank r symmetric tensors ([3]). 

The same ambient space contains the fradeco variety 71,n,d and all its secant varieties (JsTr,n,d- 

Theorem 4.7. For any r > n > d > 2, we have 

C Tr,n,d C , (30) 


and hence Tr-n,n,d C 71,n,d whenever r > 2n. Also, if r = rir 2 with ri > 2 and r 2 > n, then 


a. 


iTr2,n,d ^ 7r^n,d- 


(31) 
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Proof. We fix d. The right inclusion in fl30|) is immediate from the dehnition. For the left 
inclusion we use the parametrization of J^r,n given in f[T5D and fITB]) . The point is that the 
(r — n) X n-matrix W can be chosen freely. Equivalently, the projection of Qr,n C to 

any coordinate subspace is dominant. This means that the hrst r — n summands 

in ([1]) are arbitrary powers of linear forms, and this establishes the left inclusion in fl30|) . 

To show the inclusion ([3T]), we consider arbitrary frames Vi, V 2 , • • •, K-i ^ J^r 2 ,n- Then the 
nxr-matrix V = (Vi, V 2 ,..., Vj-i) is a frame in Each Vi together with a choice of Aj G 
determines a point on Tr 2 ,n,d- Thus we have ri points in Tr 2 ,n,d, and any point on the P’’i“^ 
spanned by these lies in Tr,n,d, where it is represented by V with A = (Ai,..., Aj-J G □ 

Example 4.8. Let n = 2 and write H = (tj+j) for a Hankel matrix of unknowns with 
r + 1 rows and sufficiently many columns. The secant variety is dehned by the ideal 

Ir+i{H) of (r+1) X (r+l)-minors of H. The ideal-theoretic version of fl30l) states that 

D Ir-l{Mr) D Ir+l{H). 

It is instructive to check this. The left inclusion follows from the Cauchy-Binet Theorem 
applied to = A ■ H where A is the (r—1) x (r-|-l) integer matrix underlying Mr. <0 

Remark 4.9. (a) Since concatenations of frames in M” are always frames, fl3T|) generalizes 

from secant varieties to joins. Namely, if r = ri -|- r 2 , then % 2 ,n,d C Tr,n,d- 

(b) The inclusion in is always strict, with one notable exception: <7272, 2 ,d = 71, 2 ,d- 

Theorem 14.71 implies that the Veronese variety is contained in the fradeco variety 

71,n,d with r > n. This is illustrated in Example 14.51 where we wrote the quadric that vanishes 
on 71,3,4 as a linear combination of the binomials that dehne C P^"^. The formula 

fl2^ shows that this cubic vanishes on cr 2 Z 23 P^. Similarly, we can verify that the cubic in 
Example 14.61 vanishes on (T 2 Z/ 4 P^ by writing it as a linear combination of the 3 x 3-minors 
Cijk,imm of the 6 X 6-catalecticant C matrix in fl3T|) . One such expression is 

50(7012,012 — 30(7o12,123 + 50(7o 12,034 — 30(7oi2,125 + 50(7 o12,045 + 63(7o12,345 — 10(7o13,024 + 10(7o13,234 
+5(7oi3,o15 + 35(7oi3,135 -|- 34(7oi3,245 + 5(7o23,023 — 80(7o23,134 + 5(7o23,025 “ 26(7o23,235 “ 19C'o23,145 
— 30(7i23,123 + 29(7i23,125 “ 10(7i23,345 “ 10(7o14,025 + 19(7o14,235 “ 53(7o14,145 “ 30(7o 24,245 + 5(7o34,034 
+26(7 o34,045 + 5(7o 34,345 + 50(7i34,i34 -|- 50(7 i 34,235 + 30(7i34,i45 -|- 30(7234,245 + 5(7oi5,oi5 + 26(7oi5,l35 
+50(7 oi5,245 — 5(7o25,235 “ 10(7o25,145 “ 10(7i 25,345 “ 4(7o35,035 + 5(7i35,i35 -|- 50(7 i35,245 + 5(7235,235 

-|-5(7o45,045 + 5(7o45,345 + 50(7245,245- 

Theorem IT71 gives lower bounds on the degrees of the equations defining fradeco varieties: 

Corollary 4 . 10 . All non-zero polynomials in the ideal ofTr,n,d must have degree at least r — n-\-l. 

Proof. The ideal of the Veronese variety contains no linear forms. It is generated by 2 x 2 

minors of catalecticants. A general result on secant varieties |23l Thm. 1.2] implies that the ideal 
of (7r-ni^dP”~^ is zero in degree < r — n. The inclusion C Tr,n,d yields the claim. □ 

In Table [2] we see that 71,3,4, 71,3,5, 71,3,4, 71,3,5 and 71,3,5 have their first minimal generators in 
the lowest possible degrees. However this is not always the case, as shown dramatically by 77 , 4 , 3 . 
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5 Numerical Recipes 

Methods from Numerical Algebraic Geometry (NAG) are useful for studying the decomposition of 
tensors into frames. Many of the results on fradeco varieties Tm^d reported in Sections [3] and 0] 
were discovered using NAG. In this section we discnss the relevant methodologies. Our experiments 
involve a mixture of using Bertini [2], ]yiacaulay2 [13], Maple, and Matlab. 

All algebraic varieties have an implicit representation, as the solution set to a system of poly¬ 
nomial equations. Some special varieties admit a parametric representation, as the (closure of the) 
image of a map whose coordinates are rational functions. Having to pass back and forth between 
these two representations is a ubiquitous task in computational algebra. 

The fradeco variety studied in this paper is given by a mixture of implicit and parametric. Our 
point of departure is the implicit representation Q of the funtf variety or its homogenization 
Gr,n- Built on top of that is the parametrization ([T]) of rank r tensors: 


(^nxr ^ 

Sym,(C-) 


U 

U 


X C” 

^ ' r,n,d 

(32) 

(^,A) 

^ Aivf+ Aavf ^ + • • • + Xrvf^ 



Here, Tr^n,d denotes the affine cone over the projective variety Tr^n,d- The inpnt to our decomposition 
problem is an arbitrary symmetric nxnx • • • xn-tensor T and a positive integer r. The task is to 
decide whether T lies in Tr,n,d, and, if yes, to compute a preimage (V, A) under the map in (l32]i . 
Any preimage must satisfy the non-trivial constraint V G 

5.1 Decomposing fradeco tensors 

We discnss three approaches to finding frame decompositions of symmetric tensors. 

5.1.1 Tensor power method 

Onr original motivation for this project came from the case r = n of odeco tensors [2T]. If T G Tn,n,di 
then the tensor power method of [Tj reliably reconstructs the decomposition ([T|) where {vi,..., v„} is 
an orthonormal basis of M”. The algorithm is to iterate the rational map VT : given 

by the gradient vector VT = {dT/dxi,... ,dT/dxn)- This map is regular when the hypersurface 
{T = 0} is smooth. The fixed points of VT are the eigenvectors of the tensor T. Their number was 
given in [7|. The punchline is this: if the multipliers Ai,..., in ([T|) are positive, then vi,..., v„ 
are precisely the robust eigenvectors, i.e. the attracting fixed points of the gradient map VT. 

This raises the question whether the tensor power method also works for fradeco tensors. The 
answer is “no” in general, but it is “yes” in some special cases. 

Example 5.1. Let n = 2, r = 4, d = 5 and consider the fradeco qnintic 

T = ax^ + y^ + {x + yf + {x - yf G 7^, 2,5, 

where a > 6 is a parameter. The eigenvectors of T are the zeros in of the binary quintic 

dT ^ ^, 4/0 1 2^2 , 31 4 

^ 5y((ax-6)x + [2xy --y ) + —y 
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The point (1 : 0) is an eigenvector, but there are no other real eigenvectors, as the expression is a 
sum of squares. Hence the frame decomposition of T cannot be recovered from its eigenvectors. <0> 

Example 5.2. For any reals Ai, A 2 , A 3 , A 4 > 0 and any integer d > 5, we consider the tensor 

T = Ai(-5a;i + X 2 + + A 2 (xi - 6 x 2 + xs)'^ + A 3 (xi + X 2 - 6 x 3 )'^ + A4(3xi + 8 x 2 + 8 x 3 )'^. ( 88 ) 

This tensor has precisely four robust eigenvectors, namely the columns of the matrix H in ([7|). 
Hence the frame decomposition of T can be recovered by the tensor power method. <0> 

The following conjecture generalizes this example. 

Conjecture 5.3. Let r = n + 1 < d and T G Tn+i,n,d with Ai,...,A,i+i > 0 m (Qp. Then 
vi,..., v„_|_i are the robust eigenvectors ofT, so they are found by the tensor power method. 

Example 15.11 shows that Conjecture 15.81 is false for r > n + 2, and it suggests that the Tensor 
Power Method will not work in general. We next discuss two alternative approaches. 

5.1.2 Catalecticant method for frames 

The matrices in Theorem 18.11 furnish a practical algorithm for the decomposition problem when 
n = 2. This is a variant of Sylvester’s Catalecticant Algorithm, and it works as follows. 

Our input is a binary form T G Tr, 2 ,n- We seek to recover the tight frame into which T 
decomposes. Since we do not know the fradeco rank r in advance, we start with AI 3 , AI 4 , AI 5 , etc. 
and plug in the coordinates U of T. The fradeco rank is the first index r with Air rank deficient. 

If the matrix Air is rank deficient, then its rank is at most r — 2. Let us assume that the rank 
equals exactly r — 2. Otherwise T is a singular point (cf. Remark [83]). Then, up to scaling, we find 
the unique row vector w G in the left kernel of Air- By Theorem 18.11 we know that Air is 

the product of the matrix Mr and an (r — 1 ) x (d — r — 1 ) matrix with entries where 

V = (vij) G Qr ,2 is the desired frame. Moreover, the matrix Mr has rank r — 2, so the vector w 
also lies in the left kernel of Mr, i.e. w • Mr = 0. Thus, 

0 = wMr = {f{vu,V2l),f{vi2,V22),---,f{vir,V2r)), 

where f{x,y) is a binary form of degree r. The r roots of f{x,y) in are the columns of the 
desired V = (vij) G Gr, 2 - Using these Vij, the given binary form has the decomposition 

r 

T{x,y) = ^Aj(uijX + U 2 j 1 /)‘^, 
i=i 

where the multipliers Ai,..., A^ are recovered by solving a linear system of equations. 

Example 5.4. Let r = 5 and d = 8. We illustrate this method for the binary octic 

T = (-287-896a)x8 + 8(65+241a)x^y - 28(16+68a)x®y2 + 56(5+8lQ;)x5y3 70(2-56a)x'^y^ 

+56(-7 + 198a)x3y5 + 28(82 - 7lQa)x‘^y^ + 8(-115 + 2671a)xy^ + (485 - 9968a)/, 

where a = y/S — 2. The parenthesized expressions are the coordinates to,... ,ts. We find 

/-18548a + 595 8686 a - 150 -996a + 42 848a + 18 \ 

_ 2092a - 94 -548a + 26 100a - 22 148a + 50 

“ -2092a+ 94 548a - 26 -100a+ 22 -148a - 50 

V 996a-80 -848a - 6 896a + 90 -1286a - 817/ 
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This matrix has rank 3 and its left kernel is the span of the vector w = (0,1,1,0). Therefore, 



T 

( Ax + 5'f'ii 
viiAi - 

^22 + 5^12 
'^12^22 ~ ^'^12^22 

A 3 + 5'^13 
^13^23 ~ 

’24 + 5’i4 
Vl4vl^ - 31)^4 

As + ^As ^ 

’’15’’25 ~ ^"^IsAs 

1 


1 


3*^11 *^21 -Ai'^21 

^^12^22 ~ ‘^12‘^22 

3*^13 *^23 ~ ’'13’'23 

3’'i4’'24 “ *'14^24 
5’'i4''24 + Ai 

^AsAs ~ ’'i5’'25 

\oJ 


V ^AiAi + Ai 

5^12^22 +"^^12 

5’'i3’'23 + ’'W 

+’'15 / 


Hence the five columns of the desired tight frame V = (vij) are the distinct zeros in of 
V2i) = Vuvti - - v\^V2i for i = 1,..., 5. 


fo = 


We hnd 

'1 0 1 a 1\ ^ 

,0 1 1 1 aj ^ 

It remains to solve the linear system of nine equations in A = (Ai,..., A 5 ) given by 
T = Xix^ + X2y^ + Xsix + y)^ + X^^ax + y)^ + X^i^x + ay)^. 

The unique solution to this system is Ai = A 2 = A 3 = A 5 = 1 and A 4 = 1552 + 896\/3. 


❖ 


5.1.3 Waring-enhanced frame decomposition 


We now examine the decomposition problem for n > 3. Since no determinantal representation 
of Tr^n^ is known, a system of equations must be solved to recover (fo A) from a given tensor in 
Tr,n,d- In some special situations, we can approach this by taking advantage of known results on 
Waring decompositions. For instance, in Example 11.11 the Waring decomposition is already the 
frame decomposition. Example 13.51 shows that this is an exceptional situation. 

We demonstrate the “Waring-enhanced” frame decomposition for the ternary quartic 

J2i+j+k=4 = 467a:^-bl52x3?/-M448x3z-h660x2y2_i4883;2yz-p4020x2z2-p536xy3 

— 1992xy^2;-|-2352xyz^-|-944xz^-|-227y^—1000y^2;-|-2148y^2;^—1960y^^-|-12672;^. 


Ternary quartics of rank < 5 form a hypersurface of degree 6 in The equation of this hyper¬ 
surface is the determinant of the 6x6 catalecticant matrix C. Here the dimension is one less than 
expected; this is the first entry in the Alexander-Hirschowitz list (|27p . Eor the given quartic. 


0 

0 

folO 

fool 

fo20 

foil 

fo02 


■ 467 

38 

362 

no 

-124 

670 ■ 

foio 

fo20 

foil 

fo30 

^121 

fol2 


38 

110 

-124 

134 

-166 

196 

fool 

foil 

^202 

fo21 

fol2 

fo03 


362 

-124 

670 

-166 

196 

236 

fo20 

fo30 

fo21 

fo40 

fo31 

t022 


no 

134 

-166 

227 

-250 

358 

foil 

fo21 

fol2 

fo31 

t022 

fol3 


-124 

-166 

196 

-250 

358 

-490 

.^202 

fol2 

fo03 

t022 

fol3 

0 

0 


_ 670 

196 

236 

358 

-490 

1267_ 


(34) 


This matrix has rank 5 and its kernel is spanned by the vector corresponding to the quadric 
q = 14m^ —uv — 2uw — 4v^ — llvw — lOic^. The points {u : v : w) in P^ that lie on the conic {q = 0} 
represent all the linear forms ux + vy + wz that may appear in a rank 5 decomposition. 

Our task is to find five points on the conic {q = 0} that form a frame V S ^ 5 , 3 . This translates 
into solving a rather challenging system of polynomial equations. One of the solutions is 

/-I 2 2 I + 2V3 -l + 2^/3\ 

V = (vi, V2,V3,V4,V5) = 2 2-1 -2 + V^ 2 + y/3 \ . 

\ 0 1 -2 5 -5 / 


The given ternary quartic has the frame decomposition -|- v®^ -|- v®^. 
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5.2 Exploring the fradeco variety 

The following tasks make sense for any variety X C arising in an applied context: (i) sample 
points on X, (ii) compute the dimension and degree of X, (iii) compute an irreducible decomposition 
of X, (iv) hnd a parametrization of X, (v) hnd some polynomials that vanish on X, (vi) determine 
polynomials that cut out X, (vii) hnd generators for the ideal of X. Numerical algebraic geometry 
(NAG) furnishes tools for addressing these points. In our study, X is the fradeco variety %^n,d- We 
used NAG to hnd answers in some cases. In what follows, we explain our computations. Particular 
emphasis is placed on the results reported in Section 0] for the degree and Hilbert function of Tr^n,d- 
All computations are carried out by working on the affine cone Tr^n,d C Symrf(C”). 

5.2.1 Dimension and degree 

The dimension and degree of the affine variety Tr^n,d can be computed directly from the mixed 
parametric-implicit representation in (I32p . The dimension can be found by selecting a random 
point on x M’', determining its tangent space via [26] . and then taking the image of this 
tangent space via the derivative of the map The image is a linear subspace in Sym^(M”’), 
and its dimension is found via the rank of its dehning matrix. These matrices are usually given 
numerically, in terms of points sampled from J>,n) so we need to use singular value decompositions. 

The computation of the degree is carried out using monodromy. We obtained the results of 
Theorem 14.31 by applying essentially the same technique as in nail], adapted to our situation 
where the mapping is from an implicitly dehned source. Here are some highlights of this method 
for Tr^n,d- We performed these computations using Bert ini and MatLab. 

Let c denote the codimension of Tr^n,dj as given by the formula in Conjecture 14.21 The degree 
of Tr,n,d is the number of points in the intersection with a random c-dimensional affine subspace of 
Sym^(C”). Here we represent the fradeco variety purely numerically, namely as the set of images 
of points (H, A) under the parametrization shown in (j32p . This method verifies the dimension 
of 7r,n,d because the intersection would be empty if the dimension were lower than expected. 

As a first step, we compute a numerical irreducible decomposition of the funtf variety This 
also gives its degree and dimension, as shown in Table [TJ In particular, we obtain degree-many 
points of J-r^n that lie in a random linear space of dimension equal to codim(J>^„). 

We take V to be one of these generic points in we select a random vector A G C^, and 
we compute the fradeco tensor A). We also fix a random c-dimensional linear subspace TZ of 

Sym^(C”) and a random point U in the c-dimensional affine space TZ + U. 

By construction, the affine cone Tr^n,d and the affine space TZ + U intersect in deg{Tr^n,d) many 
points in Sym^(C"'). One of these points is T,d{V, A). Our goal is to discover all the other intersection 
points by sequences of parameter homotopies that form monodromy loops. Geometrically, the base 
space for these monodromies is the vector space quotient Sym^pC^j/i?. 

We hx two further random points Pi and P 2 in Sym^(C"'). These represent residue classes 
modulo the linear subspace TZ. The data we hxed now dehne a (triangular) monodromy loop 


(TZ + U) r\ Pr,n,d 



(TZ + P 2 ) n Tr,n,d ^+ .Pi) G Tr,n,d 
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We use Bertini to perform each linear parameter homotopy. This constructs a path (Vt, Xt) in the 
parameter space. Here t runs from 0 to 3. We start at (Vqj-^o) = point lies 

in {TZ + Pi) n Tr,n,d for i = 1,2, and TidiV^, A 3 ) is back in (7^ + 17) n Tr,n,d- With high probability, 
XldiVs, A 3 ) / Xld{V, A) holds, and we have discovered a new point. Then we iterate the process. Let 
Sk := {Srf(17, A),..., Tidiy', A')} denote the subset of iJZ + U) r\ %^n,d that has been found after 
k steps. In the next monodromy loop we trace the paths of Sk to produce 5fc+i, the endpoints of 
monodromy loops starting from Sk- Using MatLab, we then merge the point sets to form Sk+i = 
Sk U Sk+ 1 - We repeat this process until no new points are found after 20 consecutive monodromy 
loops. The number of points in Sk is very strong numerical evidence for the degree of Tr,n,d- At 
this point, one can also use the trace test [25] with pseudowitness sets m to confirm that degree. 

5.2.2 Numerical Hilbert Function 

We wish to learn the polynomial equations that vanish on Tr,n,d- The set I of all such polynomials 
is a homogeneous prime ideal in the polynomial ring over Q whose variables are the entries 
of an indeterminate tensor T. We write this polynomial ring as 

Q[T] = 0Q[T]e ^ Sym,(Sym^(Q-)) = 0 Sym,(Symrf(Q-)). 

e >0 e >0 

The space of all polynomials of degree e in the ideal I is the subspace 

le = InQ[T]e C Q[T]e ~ Symg(Symrf(Q")). 

A natural approach is to fix some small degree e and to ask for a Q-linear basis of le- 
The dimensions of these vector spaces are organized into the Hilbert function 

N —^ N, e !->• dimQ(/e)- 

We used Bertini and Matlab to determine specific values of the Hilbert function. In some cases, 
an independent Maple computation was used to construct a basis for the Q-vector space le- 

Fix values for r, n, d- As discussed above, we can use the parametrization to produce many 
sample points T = T,d{y,X) on Tr,n,d- The condition f{T) = 0 translates into a linear equation in 
the coefficients of a given polynomial / £ Q[T]e, and le is the solution space to these equations 
as T ranges over Tr,n,d- We write these linear equations as a matrix whose number of columns is 

dim(Q[T]e) = ), and with one row per sample point T. In practice we take enough 

sample points so that le is sure to equal the kernel of that matrix. 

This procedure may be carried out in exact arithmetic over Q when sufficiently many exact 
points can be found on Pr,n- When floating point approximations are used, some care is required in 
choosing the appropriate number of points and a sufficient degree of precision. This numerical test 
can become inconclusive in high dimension due to these issues. Using floating point arithmetic and 
30,000 points of J>,n we obtained the values listed in Table [3| The blanks indicate that we did not 
find conclusive evidence for the exact value of dim(/e) in that case. For 75 , 4 , 3 , 76 , 3 , 4 , 76 , 4 , 3 , 77 , 3 , 5 , 
and 78 , 3,5 we also found no conclusive numerical evidence for equations in degrees less than 5. 

The calculation of dim(7e) is a numerical rank computation via singular value decomposition, 
so at least in principle it is possible to also extract a basis of le- However, in practice, round-off 
errors yield imprecise values for the coefficients of the basis elements of le- This makes it difficult 
to reliably determine an exact Q-basis of le by numerical methods. 
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ideal 

2 

3 

4 

5 

6 

dimX(7i,2,9)e 

0 

0 

5 

46 

235 

dimX(7I,3,4)e 

6 

127 

1093 

5986 


dimX(7I,3,5)e 

27 

651 

6370 



dimX(7i,3,4)e 

0 

1 

21 



dimX(7i,3,5)e 

0 

20 

633 



dimX(7^,3,5)e 

0 

0 

1 




Table 3: Numerical computation of the Hilbert functions of fradeco varieties 

To discover the explicit ideal generators displayed in Sections [3] and 01 we instead used exact 
arithmetic in Maple. A key step was to produce points in the funtf variety J-r,n that are defined 
over low-degree extension of Q, and to map them carefully via To accomplish this, we used 
the representation of Qr^n discussed in Section [21 In our experiments, we found that the solve 
command in Maple was able to handle dense linear systems with up to 3, 500 unknowns. 
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